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Abstract: We calculate the next to the leading order Casimir effect for a real scalar field, 
within c/> 4 theory, confined between two parallel plates in three spatial dimensions with the 
Dirichlet boundary condition. In this paper we introduce a systematic perturbation expan- 
sion in which the counterterms automatically turn out to be consistent with the boundary 
conditions. This will inevitably lead to nontrivial position dependence for physical quanti- 
ties, as a manifestation of the breaking of the translational invariance. This is in contrast 
to the usual usage of the counterterms in problems with nontrivial boundary conditions, 
which are either completely derived from the free cases or at most supplemented with the 
addition of counterterms only at the boundaries. Our results for the massive and massless 
cases are different from those reported elsewhere. Secondly, and probably less importantly, 
we use a supplementary renormalization procedure, which makes the usage of any analytic 
continuation techniques unnecessary. 
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1. Introduction 

During the last fifty years many papers have been written on the Casimir effect. In this 
paper we introduce a new approach in regards to the renormalization program. We use 
our approach in one of the simplest nontrivial example possible, i.e. a real scalar field 
confined between two parallel plates in 3+1 dimensions, with </> 4 self-interaction. As we 
shall see, our results for the next to leading order term (NLO) differs significantly from 
what exists in the literature. It is therefore suitable to start at the beginning. In 1948 
H.B.G. Casimir found a simple yet profound explanation for the retarded van der Waals 
interaction [Q]. After a short time, he and D. Polder related this effect to the change 
in the zero point energy of the quantum fields due to the presence of nontrivial boundary 
conditions [g]. This energy has since been called the Casimir energy. The zero-order energy 
in perturbation theory has been calculated for various fields (see for example ||, Q). Also 
the NLO correction, which is usually called the first-order effect, has been computed for 
various fields. For the electromagnetic field this correction is said to be due to the following 
Feynman diagram , and has been computed first by Bordag and collaborators ||. 

However, note that this correction is a two loop correction in this case and is 0(e 2 ). 
Moreover the two-loop radiative corrections for some effective field theories have been 
investigated in || 0, ||]. Next, in the case of a real massive scalar field NLO correction 
to the energy has been computed in [||, [R], [ll], [T^, [l^, |l4|, |l5|, [Tql . This correction is 
a two loop correction in this case but is O(X). Moreover, N. Graham et al used new 
approaches to this problem by utilizing the phase shift of the scattering states [17], or 



replacing the boundary conditions by an appropriate potential term 15]. However, all of 
the authors who use counterterms to calculate NLO correction to the Casimir energy use 
the free counterterms (in the space between the plates) by which we mean the ones that 
are relevant to the free cases with no nontrivial boundary conditions, and are obviously 
position independent. Only in [13| the author notes that in certain cases, counterterms can 
depend on the distance between the plates. The first use of nontrivial boundary conditions 
for the renormalization programs in problems of this sort seems to be due to Fosco and 
Svaiter IS ] . These authors use free counterterms in the space between the plates and place 
additional surface counterterms at the boundaries. Later on various authors proposed the 
use of exactly the same renormalization procedure for various physical problems ppj] . The 
first calculation for the NLO of Casimir energy for the massive scalar field using this 



renormalization program is done in ref.[21|. We should note that their results for the 



massless limit in 1+1 dimensions is infinite similar to [15], who used free counterterms 
only. Moreover, the results for the massless limit in 2 + 1 case reported in ref.[15, 21], 
depend crucially on the order in which the limits d — > 2 and m — > are taken. When the 
order of the limits are as shown, their result is infinite. When the order is opposite, they 
get finite results which, surprisingly, contains the Euler-Mascheroni number for both the 
massive and massless cases. This is in contrast to the usual point of view that this number 
should not appear in any results which reflect a physical quantity [22, 23]. We interpret 
these negative results as indications that the use of free counterterms in these problems 
might not be appropriate. It is also worth mentioning that all the papers on the analogous 
calculations of the NLO corrections to the mass of solitons, that we are aware of, use free 
counterterms (see for example [24, 25, ^]). In references [25] the authors used the mode 
number cutoff introduced by R.F. Dashen (1974) [^] to calculate the NLO Casimir energy 
due to the presence of solitons. 

In this paper, we present a systematic approach to the renormalization program for 
problems which are amenable to renormalized perturbation theory, and contain either non- 
trivial boundary conditions or nontrivial (position dependent) backgrounds, e.g. solitons, 
or both. Obviously all the n-point functions of the theory will have in general nontrivial 
position dependence in the coordinate representation. This is one of the manifestations 
of the breaking of the translational symmetry. The procedure to deduce the countert- 
erms from the n-point functions in a renormalized perturbation theory is standard and 
has been available for over half a century. Using this, as we shall show, we will inevitably 
obtain position dependent counterterms. Therefore, the radiative corrections to all the 
input parameters of the theory, including the mass, will be in general position dependent. 
Therefore, we believe the information about the nontrivial boundary conditions or position 
dependent backgrounds are carried by the full set of n-point functions, the resulting coun- 
terterms, and the renormalized parameters of the theory. Our preliminary investigations 
have revealed that the main difference between our position dependent counterterms and 
the free ones are maximal for positions which are about a Compton wavelength away from 
the walls, although it is also nontrivial at other places. Moreover, in the limit of large 
plate separation our counterterms become position independent everywhere except within 
one Compton wavelength away from the walls. In this limit, the difference between our 



counterterms and the free ones, before the transverse momentum integration, approaches 
a constant value proportional to l/2m away from the walls and 1/m at the walls. How- 
ever since these terms will multiply at least one- loop expressions which are divergent, their 
difference will have significant consequences. Here we use this procedure to compute the 
first-order radiation correction to the Casimir energy for a real scalar field in 3+1 dimen- 
sions with (j) 4 self interaction. We compute this correction for both a massive and a massless 
scalar fields and show that the massless limit of the massive case exactly corresponds to 
the massless case. 

In addition, up to now most of the papers on the Casimir effect, that we are aware of, 
use some from of analytic continuation. We share the point of view with some authors such 
as the ones in [14| that the analytic continuation techniques are not always completely 
justified physically. Moreover, like the authors of the first of the aforementioned references, 
we have found counterexamples, which we point out in this paper and elsewhere [^7], p8fl . 
The counterexamples show that analytic continuation techniques alone might not yield 
correct physical results, and sometimes even give infinite results |2{|. Therefore, we prefer 
to use a completely physical approach by enclosing the whole system in a box of volume 
V = L 3 , which eventually can go to infinity, and calculating the difference between the zero 
point energies of two different configurations. The main idea of this method is actually 



due to T.H. Boyer |3C], who used spheres instead of boxes. This we shall call the "box 
renormalization scheme" and can be used as a supplementary part of other usual regular- 
ization or renormalization programs. This box renormalization scheme, has the following 
advantages: 

1. Use of this procedure removes all of the ambiguities associated with the appearance 
of the infinities, and we use the usual prescription for removing the infinities in the 



regulated theory, as explained in Sec. |3.2| . This is all done without resorting to any 
analytic continuation schemes. 

2. In order to calculate the Casimir energy we subtract two physical configuration of 
similar nature, e.g. both confined within finite regions, and not one confined and the 
other in an unbounded region. 

3. This method can be used as a check for the cases where analytic continuation yields 
finite results, and more importantly, can be used to obtain finite results when the 
former yields infinite results. 

More importantly, we have discovered for the case of parallel plates, much to our 
surprise, that when the problem is set up correctly, that is when proper counterterms 
are used, using analytic continuation technique for NLO correction it seems impossible to 
obtain finite results in any integer space-time dimensions, and correct finite results in any 
non- integer space-time dimensions. We will make this statement more explicit later on in 
Section ||. However, this technique also gives the correct results for the leading term for 
this geometry (see for example Q ) . We should mention that some authors believe that use 
of box regularization or renormalization procedures, in which the size of the box eventually 



goes to infinities could be avoided by using appropriate boundary conditions on the fields 
at spatial infinity [31]. 

In Section 2 we calculate the leading order term for the Casimir energy in d space 
dimensional case. We do this first of all to explain more completely the physical content of 
the problem and set up our notations. Secondly this computation is just about as easy to do 
in d dimensions as is in the three dimensional case. In Section 3 we compute the first order 
radiative correction to this energy. In order to do this we first state the renormalization 
condition, and then derive an expression for the first order radiative correction for a massive 
and massless scalar field. In Section 4 we discuss the validity of the analytic continuation 
techniques relevant to our problem. In Section five we give a brief summary of our results 
and state our conclusions. 



2. The Leading Term of the Casimir Effect 

The lagrangian density for a real scalar field with </> 4 self- interaction is: 

= ~ - ^(*) 4 , (2.1) 

where mo and Ao are the bare mass and bare coupling constant, respectively. Here we 
calculate the leading term for the Casimir energy in d spatial dimensions. Obviously the 
leading term, in contrast to the higher order corrections, is independent of the form of 
the self-interaction. The Casimir energy is in general equivalent to the work done on the 
system for bringing two parallel plates from ±oo to ±a/2. As mentioned before, part of our 
renormalization procedure is to enclose the whole system in a d dimensional cubical box 
of sides L. To compute this leading term, we first compare the energies in two different 
configurations: when the plates are at ±a/2 as compare to ±6/2. We name the axis 
perpendicular to the plates the z axis. To keep the expressions symmetrical, we choose the 
coordinates so that the edges of the confining box are at ±L/2 in any direction. 
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Figure 1: The geometry of the two different configurations whose energies are to be compared. 
The labels al, etc. denote the appropriate sections in each configuration separated by the plates. 



The total zero point energy of the upper configuration in figure (||) will be called 
and of the lower one E a . In our box renormalization scheme we need to define the Casimir 
energy as follows 
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where, 



E a — E ai + 2E a2 , Ef, — E^ + 2Eh 2 . 



(2.2) 



(2.3) 



Here we choose the Dirichlet boundary condition on all of the boundaries. Then we can 
expand the field operator p in the eigenstate basis appropriate to this boundary condition, 
and its explicit second quantized form, for example in region a\ becomes 
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where, 



nir 



and n = 1, 2, . . . . 



(2.5) 



Here and k aijTl denote the momenta parallel and perpendicular to plates, respectively. 
Also and a n are creation and annihilation operators obeying the usual commutation 
relations: 



[ a ri) a n'] — ^n,vJ i [ a rt) a n'] — [ a n; a n'] — 0) 

and a|0) =0 defines the vacuum state in the presence of boundary conditions. Using the 
above equations one can easily obtain 
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(2.6) 
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where TC^ denotes the usual free Hamiltonian density, easily obtained from the Lagrangian 
density, and the superscript (0) denotes the zero (or leading) order term of this energy. 

Also k = |k |, and £ld = is the solid angle in d-dimensions. Therefore, 
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where, 



g(n) — uj ai ,n + 2w a2in — w& liTl — 2a;b 2)n . 



Now we are allowed to use the Abel-Plana summation formula, since we now expect the 
summand to satisfy the strict conditions |36| for the validity of this formula. That is, we 
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expect any reasonable renormalization program for calculating any measurable physical 
quantity to yield finite results. The Abel-Plana summation formula gives 
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where g(0) vanishes in this case due to our box renormalization. The second term in the 
bracket, using suitable changes of variables, becomes 
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where k for example in the first term denotes — treated as a continuous variable. The 

a 

above calculation shows that this term is exactly zero. Therefore, only the branch-cut term 



(the last term in eq. Q2.8D) gives nonzero contribution and the final result is 
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where K n (x) denotes the modified Bessel function of order n. Using eq. (2^) for the 
Casimir energy and noting that K n (x) is strongly damped as x goes to infinity, only the 
first term remains when the limits are taken, and the result is 
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This is the result for the leading term for the Casimir energy in d-dimensions, on which all 
authors agree (see for example ||, H|). It is important to note that we, unlike most other 
authors, derived this result without any use of analytic continuation. If we set d = 3, we 
have 
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The results are in agreement with what exists in literature (see for instance |l4|]). It is 
interesting to note that for the massless case, the result is, not surprisingly, exactly half of 
the corresponding expression for the electromagnetic case. 



3. First-Order Radiative Correction 



Now we calculate the next to the leading order (two loop quantum correction) shift of the 
Casimir energy for a scalar field in (f) 4 theory using the renormalized perturbation theory 
in 3 + 1 dimensions. As mentioned before, the main idea of our work is that when a 
systematic treatment of the renormalization program is done, the counterterms needed to 
retain the renormalization conditions, automatically turn out to be position dependent. 
This, as we shall see, will have profound consequences. However, our main scheme of 
canceling the divergences using counterterms and a few input experimental parameters, is 
in complete conformity with the standard renormalization approach. To set the stage for 
the calculations, we shall very briefly state the renormalization procedure and conditions. 

3.1 Renormalization Conditions 

The <^> 4 Lagrangian eq. (|2.1| ), after rescaling the field ip = Z 1 ^ 2 Lp T) where Z is called the 
field strength renormalization, and the standard procedure for setting up the renormalized 
perturbation theory, becomes (see for example ]23]1 ), 



+\^z[d^ T {x)] 2 - l-5 m (p r (x) 2 - ^v? r (x) 4 , (3.1) 



where 5 m ,5\,5z are the counterterms, and m and A are the physical mass and physical 
coupling constant, respectively. In this problem we are to impose boundary conditions 
on the field at the walls. An alternative approach would be to add appropriate external 
potentials to the Lagrangian so as to maintain the boundary conditions on the fields |18|| . 
We will use the first approach. Obviously the presence of nontrivial boundary conditions 
breaks the translational invariance and hence momenta will no longer be good quantum 
numbers. Therefore we find it easier to impose the renormalization conditions in the 
configuration space. For example, the standard expression for the two-point function is, 
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Since the birth of quantum field theory, as far as we know, the assertion has always been 
that the above expressions can be expanded systematically when the problem is amenable 
to perturbation theory. For example, in the context of renomalized perturbation theory, as 
indicated in eq.(|3.1|), we can symbolically represent the first few terms of the perturbation 
expansion of eq. (|3.2|) by 



Q 



~0^=. . +_v-* + — e— +.... (3.3) 

Xl V_X X 2 Xi X 2 Xi X X 2 Xl X x 2 

where • — © — • refers to the appropriate counterterm. It is obvious that the above ex- 

Xl X x 2 

pression represents a systematic perturbation expansion, and most importantly, all of the 



propagators on the right hand side should be the one appropriate to the problem under 
consideration, that is they should have the same overall functional form as the first term. 
Our first renormalization condition is that the renormalized mass m should be the pole of 
the propagator represented by the first term in (3.3). This implies the second and third 
diagrams should cancel each other out in the lowest order, and this in turn implies the 
cancelation of the UV divergences in that order, and that the counterterms will in general 
turn out to be position dependent. The renormalized mass m will then naturally turn 
out to be position dependent as well. However, we only need to fix the value of m(x) at 
one position between the plates by our renormalization condition. The exact functional 
dependence of m(x) will then be completely determined by the theory. That is, we insist 
the overall structure of the renormalization conditions such as above, and the counterterms 
appearing in them should be determined solely from within the theory, and not for example 
be imported from the free case. The equations are self deterministic and there is no need 
to take such actions. Obviously we still need a few experimental input parameters for the 
complete renormalization program, such as m{x) for some x. Analogous expression and 
reasonings could be easily stated for the four-point function. 

To one-loop order, the standard renormalization conditions applied to eq. (|3.3|) and its 
four-point counterpart, give 



5 z (x) = 0, S m (x) 
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(3.4) 



respectively. Here G(x,x') is the propagator of the real scalar field and x = (t, x). Obvi- 
ously the counterterms automatically incorporate the boundary conditions and are position 
dependent, due to the dependence of the two and four-point functions on such quanti- 
ties. Now, the higher order contributions to the vacuum energy in the interval al (i.e. 
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where^g^= —i5\(x) and — — = i[p 2 6z(x) — 5 m (x)] refer to the counterterms. Accord- 
ingly, the 0(A) contribution to the vacuum energy is 
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where G ai (x,x') is the propagator of the real scalar field in region al. Using eqs.([3.4|) and 
( |3.6| ), we obtain, 
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3.2 The Massive Case 



As mentioned before, here we choose the Dirichlet boundary condition on the plates. Then, 
after the usual wick rotation, the expression for the Green's function in the four dimensional 
Euclidean space becomes 
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It is very important to note that the 5 m (x) = -—G(x,x) (c.f. eq. (p^)) is explicitly 
position dependent. Using eq. (3.8) and eq. (3.7) and carrying out the spatial integration 
one obtains, 
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To compute the integrals we regulate them by a momentum cutoff A. Then we have 
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where io' al n 2 = k 2 ^ n + m 2 . Now the NLO for the Casimir energy, eq^^), has four terms 
the first one of which is shown above. It is clear that we can adjust the cutoffs for each 
region so that the quadratic divergences coming from different regions cancel each other. 
Similar cancelation occurs for the linear divergences. The cancelation of these divergent 
quantities without any residual finite terms is the usual practice in regulated theories [ 14 ] , 
and this is the prescription that we shall use. We shall comment on this further when we 
use it one more time for eq.( p,12| ) below. Hence we get 
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This computation is obviously complicated and plagued with a multitude of infinities. As 
explained before using the usual renormalization programs in conjunction with our box 
renormalization scheme, should eliminate all of the infinities, as might be apparent from 
the above equation. However, proper regularization schemes should still be implemented 
and proper care taken when handling these infinite expressions. All the summations ap- 
pearing in the eq. ( |3.11| ) are separately infinite. We want to use the Abel-Plana formula to 
convert these sums into integrals. However these sums do not satisfy the stringent require- 
ments stated in the Abel-Plana theorem for such a conversion [37|. However our box 
renormalization scheme provides a solution by subtracting these double sums as indicated 
in eq. ([2.2|). Now we can expect this new summand to satisfy the requirements for the 
Abel-Plana theorem. Then all the infinities actually cancel and the result for the two-loop 
correction reduces to (see Appendix for more details): 
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where f(a) = B(a) ( — — J and B(a), defined by the following expression 
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refers to the so called branch-cut term of the Abel-Plana summation formula and is a 
finite quantity. Note that the last integral in eq. ( |3.12| ) seems to diverge so it must be 
properly regularized. We prefer to use a regularization scheme for this integral term which 
is analogous to the zeta function regularization for the sums. That is, we set the power of 
the integrand to ^ — 1 and in the final stage we let d approach three. Henceforth we shall 
also refer to this as zeta function regularization. We should note that this regularization is 
not sufficient to avoid the use of analytic continuation. However, as we shall see, by adding 
an auxiliary cutoff regulator we can ultimately avoid resorting to any analytic continuation. 
We thus obtain 
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The asymptotic behavior of the hyper geometric function 2-F1 Q, 1 — f , §, —K 2 ^ multiplied 
by K for the large K is, 
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Referring to eq. ( 3.12j ) and adhering to the usual prescription described for eq. ( 3.1CQ , we 
can cancel the quadratic and logarithmic divergences using our full freedom to choose four 
different cutoffs in the four distinct integration regions corresponding to {al, a2, bl, b2} 
. Hence we get, 
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Using eqs. (§J), (|1|), (|3l3|) and (|1|), we obtain 
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This is the two-loop radiative correction to the Casimir energy. This result is obviously 
finite, and we believe it could not have been obtained with any regularization or analytic 



continuation schemes in common use. Our result defers from 15, 21]. 

Two particular limits are interesting to calculate; The large mass limit, ma 3> 1, and 



small mass limit, m — > 0. In these limits eq.(3.17) becomes 
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Our massless limit differs from the analogous results that can be extracted from refs. 



1 15, gl], 32, Q by a minus sign. Figure || illustrates the graphs for the leading and the 
NLO terms for the Casimir energy as a function of a, in the massive case and its massless 
limit. 



3.3 The Massless Case 

In the massless case it is sufficient to set the pole of the propagator to zero, i.e. one can 



set m = in the eq. ( |3.8|) , hence in the eq.( |3.11 ) we have, for example u/ a n 



k 



1 One may argue that ambiguities always exist in problems where one has to subtract infinite quantities, 
and the Casimir problems certainly fall into this category. Two methods are in common use: First is the 
analytic continuation techniques which, although usually yield correct results, do not have a very solid 
physical justification and also sometimes yield infinite results. We have devoted Section 4 to this subject. 
Second is the regularization schemes, which is what we have used. In the latter category when the problem 
is regularized, one can make a systematic expansion of the quantities in question in terms of the regulators. 
Then the terms which tend to infinity when the regulators are removed and the finite terms naturally appear 
separately. See for example eqs. (3.1C , 3T5[ ). What is almost invariably done is to adjust the regulators so 
that the singular terms exactly cancel each other, i.e. without extracting any extra finite piece from the 
difference between the infinite quantities (see for example [jl^]). This is also apparent in the leading term 
for the Casimir energy in eq. ( |2.9| ) where, as explained in the Appendix, The four changes of variables are 
equivalent to choosing four different cutoffs. One could have adjusted them so that as usual the infinities 
cancel, but any finite term would remain. However, the well known answer is obtained only when there 
is no remaining extra finite term in this subtraction scheme. This is the prescription that we have used. 
However, we do believe that this is a subject that needs further study. 
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Figure 2: The leading terms for the Casimir energy and its NLO corrections as a function of a 
for the massive (m = 1) and massless cases for A = 0.1. The correction term for the massless 
case is always negative. However the correction term for the massive case has a peak of magnitude 
~ +1.57 x 10~ 9 at a ~ 3.04 and asymptotes to zero from above. 



Accordingly this equation becomes 
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We use Abel-Plana summation formula to transform the sums into integrals, then one can 
remove the divergent integrals by changing variables as explained before. Then the result 
is 
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As before the eq.( |2.2[ ) gives the Casimir energy, 
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This result is in exact agreement with the small mass limit calculated in previous subsection. 

Note that we are explicitly assuming that 5 m= o(x) / 0. we like to stress that this 
quantity should not in general be a priori set to zero. This is in contrast to the view ex- 

This is yet another important counterexample 



pressed in for example refs.[32, 33, 34, I 
for the validity of analytic continuation: As is well known the massless limit of the ana- 
lytic continuation of the mass counterterm in i^> 4 theory is zero for space-time dimensions 
bigger than two. However one cannot renormalize the massless theory without the mass 
counterterm (see for example f23fl). 



4. On the Analytic Continuation Techniques 



In the previous section we saw that enclosing the whole system in a box is very useful 
for removing all of the infinities which appear in the calculations. In this section, for 
comparison purposes, we repeat some of those steps using analytic continuation techniques. 
We can then show or point out three examples and two counterexamples for the validity 
of the analytic continuation techniques. 

The first example for the validity of these techniques is the derivation of the leading 
term mentioned in the Section |2[ The second example is the following: Starting with 
eq. ( |3.8| ) which is obviously infinite, one can regularize it by dimensional regularization: 

2r(l - d/2) ^ 
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(4.1) 



This equality holds only for d < 1, however it can be analytically continued to d > 1. 
Putting this Green's function in eq. (|3.7| ) and performing the space integration and repeat- 
ing analogous steps for the other regions, one easily obtains: 
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Putting d = 3 in eq. (4.2), eq. (3.11) is obtained. Consequently we have shown that the 
analytic continuation for the Green's function, in this case, gives the correct result. The 
third example for the validity of the analytic continuation technique is that with its mere 
use one can obtain eq. ( |3.20D from eq. ( ]3 . 1 9| ) . 

The first counterexample can be illustrated in the continuation of derivation of the NLO 



Casimir energy from eq. (4.2). Using the Abel-Plana formula we encounter the following 
term 



(1 + K 2 



^dn 



2 > 



2r(i - 1) 



(4.3) 



This term is infinite in 



and this leads to the a term which contains T(l — d/2)T( 
any integer space dimension. This is an obvious counterexample for the use of this analytic 
continuation technique. As for the second counterexample, it concerns the massless limit of 
the analytic continuation of the mass counterterm, as explained at the end of the previous 
subsection. 

5. Conclusions 



o 
o 

[\J 
o 
o 

o 
o 
o 



We have introduced a new concept in this paper. We have insisted that the renormalization 
program should completely and self-consistently take into account the boundary conditions 
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or any possible nontrivial backgrounds which break the translational invariance of the 
system. We have shown that the problem is self-contained and the above program is 
accomplishable. To be more specific, there should be no need to import counterterms from 
the free theory, or even supplementing them with the ad hoc attachment of extra surface 
terms, to remedy the divergences inherent in this theory. In general this breaking of the 
translational invariance is reflected in the nontrivial position dependence of all the n-point 
functions. As we have shown, this has profound consequences. For example in the case of 
renormalized perturbation theory, the counterterms and hence the radiative corrections to 
parameters of the theory, i.e. m and A, automatically turn out to be position dependent. In 
this regard we disagree with the authors who use the free counterterms (see the Introduction 
for actual references). Obviously we still need a few experimental input parameters for the 
complete renormalization program, such as m{x) for some x. However, the interesting 
point is that the theory then completely determines m{x). 

Secondly we have used a supplementary renormalization scheme to be used along 
side the usual renormalization program. In computations of these sorts, there usually 
appears infinities which can sometimes be removed by the usual renormalization programs 
that often contain some sort of analytic continuation. These procedures are sometimes 
ambiguous. Our scheme is simply to confine the whole physical system in a box, and 
to compute the difference between the values of the physical quantity in question in two 
different configurations. Use of this procedure removes all of the ambiguities associated 
with the appearance of the infinities, and we use the usual prescription for removing the 
infinities in the regulated theory. Using our method, we have computed the zero and first 
order radiative correction to the Casimir energy for the massive and massless real scalar 
field in 3+1 dimensions. For the zero order, our results are identical with what exists in the 
literature. However, our first order results are markedly different from those reported in 
refs.[15, 21]. Our results for the massive and massless cases are different from theirs due to 
the aforementioned conceptual differences. As we have shown our results for the massless 
case and the massless limit of the massive case are identical. However, the massless limit of 
their massive results agrees with the "exact" results obtained in [ 32 , , who set 5 m equal 
to zero in their massless cases, or their equivalent. This is our second main difference in 
approach to the problem. As mentioned before, we believe that 5 m should not be arbitrarily 
set to zero even in the massless case, since in that case the renormalization conditions can 
no longer be fully implemented, although the theory is still in principle renormalizable. 
In this regards, we like to emphasize that their results for massless case is infinite in 1+1 
dimensions. 
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A. Use of the Abel-Plana summation formula in our subtraction scheme 



In this appendix we present the details of the calculations leading to eq. ( |3.12| ). The 
Abel-Plana summation formula (see for example pTjj) is: 
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In order to obtain first-order radiative correction, just as we discussed in Section for 



leading term, we need to compute — . We rewrite the eq. (|3.11|) as, 
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where S(a,n) = (m 2 + n Ji J . Using the Abel-Plana formula eq. ( A.l ), we obtain 
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o 
o 

N) 
O 
O 

o 
o 
o 



Using eq. ( |A.1| ) again, and making appropriate changes of variables to make the integrals 
dimensionless, all the actual infinities cancel and we finally obtain, 
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It is important to note that all these cancelations are easily accomplished using our 
supplementary box renormalization scheme. On a minor note, it is interesting to note 
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that the changes of the variables leading to the cancelation of infinities are, surprisingly, 
equivalent to setting different cutoff regularizations on the upper limits of the integrals. 
Equation ( |A.4j ) is our main equation for the NLO Casimir energy, and appears in the text 
as eq. ( |3.12| ), and is analyzed further there. 
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